Acoustic Scattering Predictions for 

Complex Sources 
and Aircraft Configurations 


Douglas M. Nark 


NASA Langley Research Center 
Spring Acoustics Technical Working Group Meeting 

Hampton, VA 


April 21-22, 2015 


Outline 



• Background 

• MIT D8. 5 scattering/shielding 

- FSC Suppression Tables 

- Next Steps 

• Incident Source Directivity Effects 

- F31A31 9x15 test predictions 
— Initial TD-Fast predictions 

- Additional geometric complexity 

• Concluding Remarks 


Spring Acoustics Technical Working Group, October 21-22, 2015 



Background 




Incorporate higher-fidelity acoustic scattering/shielding predictions 
into system noise assessment of integrated configurations 
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MIT D8.5 System Noise Assessment 



• Shielding can constitute a large 
portion of Stage 4 margin 

• Work from initial D8.5 (tri-jet) 
assessment by Berton 

• ANOPP/WING module initially employed 

• Create ANOPP suppression table using FSC 

• Quantify effect on system noise benefit 
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MIT D8.5 System Noise Assessment 

Modeling Progression 




FSC Modeling: 

• Incorporates geometric effects 

• Propulsion sources considered 


separate point sources 
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MIT D8.5 System Noise Assessment 

Suppression Table Generation using FSC 

AdB 

• Suppression Factor, 5” = 10 10 

• AdB ^P^shielded ^^^unshielded 

• S < 1 indicates suppression 

• S > 1 indicates amplification 

• Predictions combined for each 1/3-octave band 
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MIT D8.5 System Noise Assessment 

FSC Suppression Tables 


/= 500 Hz 


Amplification 


/ = 2000 Hz 



S: 0 0.25 0.5 0.75 1 1.25 1.5 1.75 2 2.25 2.5 2.75 3 
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MIT D8.5 System Noise Assessment 

Total Acoustic Field:/ =500 Hz, M = 0.23 
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D8.5 Assessment: Next Steps 

System noise assessment using FSC suppression tables 

• Improved estimate of shielding effectiveness 

Couple with work on adaptive low-drag liners 

• Assess system benefit with updated suppression tables (accounting for liners) 

• Adaptation on impedance and/or drag metrics 

• Optimize liner location based on source and liner characteristics 
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Incident Source Directivity Effects 

• Consider incident acoustic field from various propulsion sources 

• 14x22 HWB test: Broadband Engine Simulator (BENS) 

• Podded engine configuration: Turbofan source 

• Open Rotor: F31A31 Historical Baseline 





Van Zante et. ol ISABE-2011-1310, Stephens and Envia: AIAA Paper 2011-2940, Envia: CMFF '12 
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F31A31 9x15 Test 


Blade Loading 

■ OVERFLOW, FUN3D 

■ RPM: 6436 

■ Time Resolution: 1-degree 

Acoustics 

. ASSPIN, PSU-WOPWOP. F1A 
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Isolated 



Microphone 

Locations 

ii 


A 


F31A31 9x15 Test 



TD-Fast (Initial Set-Up): Short Barrier, Point Source 
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F31A31 9x15 Test 


TD-Fast (initial Set-Up): Long Barrier, Point Source 

{f= 1600 Hz) 


Forward Position 


Aft Position 



F31A31 9x15 Test 


Additional Complexity 
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F31A31 9x15 Test 



Full Test Section 


• Incorporate full tunnel geometry 
in scattering predictions 



Source Surfaces 
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Concluding Remarks 



• Apply higher-fidelity scattering methods to quantify effects 
on system noise assessments 

• Investigate benefits of properly placed low-drag external 
liners on system noise 

— Broadband attenuation capabilities increase possible liner locations 

• Incorporate tunnel geometry in scattering predictions 


Incorporate higher-fidelity acoustic 
scattering/shielding predictions into system noise 
assessment of integrated configurations 
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Background 


Conventional Tube and Wing 
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Suppression Table Creation 



supp: 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 

Installation Effects 
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Scattering/Shielding Methods 

ANOPP: WING module 

• Zorumski, W., "Aircraft Noise Prediction Program Theoretical Manual/' NASA TM-83199, 1982. 

DIM3 (based on MIT diffraction integral method) (era: pi: Spakovszky, poc: Burley) 

• Colas, D., "A Diffraction Integral Based Turbomachinery Noise Shielding Method/' PhD 
dissertation, Massachusetts Institute of Technology, 2011. 

Fast Scattering Code (FSC) (SFW NRA (ended 2012): PI: Dunn, POC: Nark) 

• Dunn, M., and Tinetti, A. F., "Aeroacoustic Scattering Via the Equivalent Source Method," AIAA 
2004-2937. 

• Tinetti, A., and Dunn, M., "The Fast Scattering Code (FSC): Validation Studies and Program 
Guidelines," NASA CR-2011-217158. 

Time Domain Fast Acoustic Scattering Toolkit (TD-FAST) ( as nra: pi: Fang Hu, poc: Nark) 

• Hu, F., "An efficient Solution of Time Domain Boundary Integral Equations for Acoustic 
Scattering and its Acceleration by Graphics Processing Units," AIAA paper 2013-2018. 

• Hu, F., "Further Development of a Time Domain Boundary Integral Equation Method for 
Aeroacoustic Scattering Solutions," AIAA paper 2014-3194. 
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Diffraction Integral Method (DIM3) 

• Based on Kirchoff diffraction theory expressing diffracted field as 
superposition of waves emitted through an aperture 

• Determine outline of shielding object based on source location 
and numerically solve the contour integral 

. Babinet's Principle with wedge potential for edge modeling 
. Beam tracing for reflections 
. Allowance for directional sources 
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Fast Scattering Code (FSC) 

• Frequency domain solve of a 3-D Helmholtz boundary value problem via 
the equivalent source method (ESM) 

• Scattered acoustic pressure field expanded into a series of point sources 
(Ns) distributed on a fictitious surface placed inside the actual scattering 
surface 



• Scattering surface is discretized into collocation points (Nc) to produce a 
dense, over-determined system of linear equations of size Nc x Ns. 

• Source strengths are adjusted to satisfy surface boundary condition 
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Time Domain BEM (TD-FAST) 




• Time domain boundary integral equation (TDBIE) reformulated for the 
convective wave equation 

• Numerical instability in time marching stages addressed via Burton- 
Miller type formulation 

• Computational cost of direct solution addressed via 

- High-order basis functions 

- Multi-level Fast Methods 

• Utilize multi-core CPU and GPU architectures 
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Validation Studies 




* Source projects geometrically similar shadow zone 

* Spheroids (centered at origin) 

• Sphere: r = a = b = 5.0 m 

• OS1: a = 5.77 m, b = 1.147 m 

• OS2: a = 5.77 m , b = 0.38 m 

* Sound Source 


Point source 



30 m 


• Monopole of unit strength 

• Frequencies: 1 < ka < 400 

• Observer fields 

• Bisecting plane, plane at z = -30 m 

• Line at z = -30 m 

• Ring, r = 7.5 m, centered at origin 
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Validation Studies 



Sphere/Spheroid Predictions 


DIM 


FSC 




TD-Fast 



Pope, Burley 


Tinetti (NRA: NNL09AA17C) 


Hu (NRA: NNX11AI63A) 
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A 


Validation Studies 




* Configuration geometry/data from 

* Cylinder (centered at origin) 

• Diameter: d = 0.48 m 

• Length: L = 3.05 m 

* Flat plate 

• Square edges 

• W= 0.5 m, L = 1.6 m, t = 0.07m 

* Sound Source 

• Monopole of unit strength 

• Location: (0.0, 0.0, 0.9936) m 

• Excitation frequencies: 9 < kd < 69 

* Observer fields 

• Bisecting plane 

• Line at z = -5.04 m 

• Sphere: r = 2.5 m, centered at origin 
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Validation Studies 



Fast Scattering Code 
(f = 8kHz) 

Cylinder Alone Cylinder- Plate 



Tinetti (NRA: NNL09AA17C) 
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Validation Studies 


Time Domain BEM: TD-Fast 

(Plate Alone: f = 8kHz) 
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Hu (NRA: NNX11AI63A) 
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Isolated OR Noise Prediction 
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Isolated OR Noise Prediction 
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Isolated OR Noise Prediction 



CQ 



100 

cuo 


_c 

90 

-3 


05 

O 

80 

J 

70 


60 


50 


40 


30 


o 


10 



- PSW (10 Rev) 

o 

ASP (F31) 

o 

ASP (A31) 

V 

LINPROP 



20 30 

Shaft Order 


40 


50 


Spring Acoustics Technical Working Group, April 15-16, 2014 


Isolated OR Noise Prediction 
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